Kondo Stripes in an Anderson-Heisenberg Model of Heavy Fermion Systems 
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We study the interplay between the spin-liquid and Kondo physics, as related to the non-magnetic 
part of the phase diagram of heavy fermion materials. Within the unrestricted mean-field treatment 
of the infinite-C/ 2D Anderson-Heisenberg model, we find that there are two topologically distinct 
non-degenerate uniform heavy Fermi liquid states that may form as a consequence of the Kondo 
coupling between spinous and conduction electrons. For certain carrier concentrations the uniform 
Fermi liquid becomes unstable with respect to formation of a new kind of anharmonic "Kondo 
stripe" state with inhomogeneous Kondo screening strength and the charge density modulation. 
These feature are experimentally measurable, and thus may help to establish the relevance of the 
spin-liquid correlations to heavy fermion materials. 
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The intermetallic heavy fermion compounds based on 
either rare earth elements or on actinides are prototypi- 
cal strongly correlated systems [l|, y, |3[ . In these materi- 
als, there are two types of electrons: dclocalized con- 
duction electrons (c-electrons), which derive from the 
outer atomic orbitals, and strongly localized /-electrons 
that singly occupy the inner orbitals. The interplay 
between the c-f hybridization and the screened on-site 
Coulomb repulsion, responsible for the single occupancy 
of / orbitals, gives rise to wide range of behaviors, in- 
cludingmagnctic ordering, Fermi liquid and non- Fermi 
liquid y. The magnetism arises primarily through the 
spin-ordering in the localized /-band. The heavy Fermi 
liquid (HFL) in the non-magnetic phase is believed to oc- 
cur via Kondo hybridization of the / and c bands. The 
non-Fermi liquid behavior emerges in the vicinity of the 
quantum critical point separating these phases. 

Recently, Senthil et al. [5|, l6( proposed a new route to 
a HFL state. They considered the possibility that due to 
the frustrated magnetic interaction between spins in the 
/ band, instead of forming an ordered magnetic state at 
low temperatures, a magnetically featureless spin liquid 
emerges with spin-1/2 chargelcss fermionic excitations 
(spinons). Upon inclusion of strong enough Kondo in- 
teraction between the spinons and conduction electrons, 
a HFL would form, with spinons acquiring an electri- 
cal charge [7[ and merging with conduction electrons to 
form a "large" Fermi surface. An important distinction 
between this scenario and the more conventional one [8| , 
which does not include the possibility of a spin liquid, is 
that the transition occurs at a finite value of the Kondo 
coupling. This results in a new kind of scaling behav- 
ior near the quantum critical point separating the FL* 
phase of two decoupled - electron and spinon - Fermi 
liquids and the HFL phase of coupled spinons and elec- 
trons p, lZI, |9| . Whether such a quantum phase transition 
is realized in nature, or always preempted by a conven- 
tional Landau phase transition [lO[, remains to be seen. 

Spin liquid phases have proven to be very elusive as 



the ground states of real materials, as well as microscopic 
models, particularly ones involving real electrons. How- 
ever, at finite temperatures, in systems with frustrated 
magnetic interactions, magnetically disordered phases 
are well approximated by spin liquids, even if at very low 
temperatures a small magnetic order parameter emerges. 
Therefore, it is natural to ask if the spin liquid tenden- 
cies can manifest, if not in the zero-temperature phase 
transition itself, in the structure of the "ordered" state 
away from the phase transition. In this Letter, we an- 
alyze in detail the structure of the heavy Fermi liquid 
that emerges out of the coupling of the spin liquid and 
the conduction Fermi lic^uid in the Anderson-Heisenberg 
lattice model. Our main findings are: First, that start- 
ing from the same spin liquid state one can obtain more 
than one kind of uniform heavy Fermi liquid. Second, 
that the spinon-electron Kondo instability can be gener- 
ically inhomogeneous. The resulting structure has inho- 
mogeneous distributions of Kondo couplings, as well as 
anharmonic charge density modulation, very similar to 
the "stripes" in the Hubbard model Uj. A harmonic 
analog of this phase (without charge modulation) was 
recently found in a continuum model by Paul et al. [12l |. 
We start from the Anderson-Heisenberg [13| lattice 
model with localized /-band. 
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Here Cic {c\^) annihilates (creates) a conduction electron 
with spin a on site i; fia ifi„) represents the /-electron. 
The number operators for c and / orbitals with spin a 
are given by n,^ = cl^da and h{^ = fLfi'^^ respectively. 
The quantity i^ is the hopping integral of the conduc- 
tion electrons, and e/ is the local energy level of the / 



orbitals. The hybridization between the conduction and 
/ bands is represented by Vij . The /-electrons experience 
the Coulomb repulsion of strength U. We have included 
explicitly the superexchange (Heisenber) interaction be- 
tween /-electron spins, S^ = \Y.ap fla'^aphp with a 
the Pauli matrix. The chemical potential /z is introduced 
to tune the mismatch between the c and / patches of 
Fermi surface. 

For simplicity, we study here the [/ — > oo limit, which 
excludes the double occupancy of the /-sites. This limit 
explicitly breaks the particle-hole symmetry in the / 
band. We find however, that results do not qualita- 
tively change for the Kondo-Heisenberg model [6| , which 
preserves the symmetry. The constraint can be treated 
by means of the "slave-particle" representation [ij, [l5| 
fiij ~ fiij^\ which replaces the original fermion with two 
auxiliary particles - a fermion / (spinon) and a boson 
h (holon), satisfying the local constraint: ^^fj^fia + 
6j&i ~ 1. Further, decoupling the superexchange inter- 
action via the Hubbard-Stratonovich transformation, the 
Hamiltonian ([T]) can be written as 
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The single occupancy constraint on each site has been 
implemented through a Lagrange multiplier A^. To pro- 
ceed, we make a static approximation by assuming that 
the slave bosons are frozen (bi —^ (bi)) and the spin liq- 
uid parameters assume their mean field values, Xij = 
{JH/'2){fJ^fjcr)- Moreover, for the purpose of present 
discussion we assume that the uniform spin liquid (xij ^ 
X uM) is a good reference state for the formation of the 
Kondo Fermi liquid, thereby making x Stii input param- 
eter of our mean field model. Note that there is a local 
C/(l) gauge freedom present in this description: A change 
of phase of the slave boson on site i, bi -^ biC^"^' , with the 
simultaneous change ft — * /ie**^' and Xij ~^ e~^'^^Xij^'^^ 
leaves the physical state unaltered. There are alternative 
mean-fields that are possible, including anomalous decou- 
pling of the superexchange term, similar to the t-J model 
for cuprate superconductors [l7|. However, since our fo- 
cus here is on possible inhomogeneous Kondo phases, we 
defer this possibility to a future consideration [18| . In the 
lattice space, the Bogliubov de-Gennes (BdG) mean-field 
equations are: 
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subject to the constraints 
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Here h^^ = -tlj - fiS^j, X, ,.j„^, ..^ 

—Xij + (e/ + Ai — ^)5ij. We solve this set of equa- 
tions self-consistently via exact diagonalization on a two- 
dimensional square lattice. For simplicity, we consider 
on-site hybridization only, Vij = VSij . We note, that 
a finite-ranged hybridization can lead to non-zero or- 
bital momentum Kondo states with the hybridization gap 
anisotropic in momentum space [l9[. The momentum- 
space inhomogencity is complementary to the real space 
one discussed here. Throughout this work, the quasipar- 
ticle energy is measured with respect to the Fermi energy 
and the energy unit t'^ = 1 is chosen. 




FIG. 1: (Color) A cartoon of possible Fermi surface topologies 
of conduction electrons and spinons. (a) The effective mass 
of spinons is infinite [x — 0) and no spinon Fermi surface 
is present; (b) X 7^ 0: The Fermi surfaces of the conduction 
electrons and spinons, in general, are mismatched. 

Within the mean-field treatment outlined above, there 
are two distinct uniform heavy Fermi liquid states. Sup- 
pose that we fix the gauge such that the spin-liquid pa- 
rameter is positive, x > 0, and solve for Kondo hybridiza- 
tion, bi. One possible solution is a constant, bi ~ b^. 
Another is bi = 6o(—l)*"= "*"*", which naively appears to 
break translational invariance. However, upon a gauge 
transformation it can be transformed to a uniform state 
with bi = &o and x < 0, which is clearly uniform. The 
two states can be distinguished based on sign{biXijbj) 
[20]. We will denote the Kondo phase with hxijbj > 
as "even" (eHFL) and hxijbj < as "odd" (oHFL). For 
these two states, we can solve the above set of equations 
in momentum space. For eHFL, the energy dispersion 
has two branches given by: 



Ei,,±^Zi,±JQ. 
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where Zk(Qk) = (^k ± ^k)/2 with ^^ = -2t=(cosfc^ + 
cosky) — fi and ^^ — ~2x{coskx -f-cosfcy) + (A — /i + e/). 




FIG. 2: (Color) Dispersionless spinons, x ~ 0. (a) The slave 
boson 6-field and the effective /-site energy e/ -I- A — /i, as a 
function of V for xo = M = 0. (b) The projected conduction 
(black line) and f-electron (red line) density of states, as a 
function of energy for V = 1.0. The temperature T = 0.01. 



and Ak = Vbo. For oHLF, one only needs to replace 

X^ -X- 

In Fig. [l] we provide a cartoon of possible Fermi 
topologies of decoupled conduction electrons and spinons. 
When the spinon energy is dispersionless, the Fermi sur- 
face is undefined, i.e. all spinons have the same - zero ~ 
energy in the Brillouin zone, Fig. [TJa). This is the situ- 
ation of the standard Anderson lattice model. When the 
spinons disperse in the presence of a uniform spin liquid, 
the spinon Fermi surface is formed, Fig.llfb). Again, due 
to the single occupancy constraint, it is pinned at half- 
filling, regardless of the chemical potential. The Fermi 
surface of the conduction electrons, FSc, can be tuned 
by the chemical potential in both cases. 

First, as a point of reference, we consider the case of 
dispersionless spinons, % = 0. In Fig. [DJa) we show the 
slave boson field (effective c-/ band hybridization) and 
the mean-field /-site energy, as a function of the coupling 
parameter V . In Fig.[2l^b), we show the projected density 
of states (DOS). There is a pronounced hybridization gap 
that appears above the Fermi energy in both projected 
conduction and /-electron DOS. 

Now we turn to the case of dispersive spinons. For 
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FIG. 3: (Color) The slave boson 6-field and the effective /-site 
energy ef + X — j^i for both the even and odd Kondo phases, as a 
function of V for /i = but x = 0.1. Inset: The corresponding 
Helmholtz free energy. The temperature T = 0.01. In this 
case, the odd Kondo state has a lower energy than the even 
Kondo state. 
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FIG. 4: (Color) The projected conduction (black lines) and 
/-electron (red lines) densities of states for the even (a) and 
odd Kondo (b) states, corresponding to Fig. [3] 



zero chemical potential /i = 0, the uncoupled Fermi sur- 
faces of the conduction electrons, FSc, and spinons, FS/, 
coincide, i.e., the red line in Fig. [TJb) overlaps the blue 
one. In this case, there is a perfect nesting across the two 
Fermi surfaces. This has important consequences for the 
Kondo hybridization, breaking the degeneracy between 
the odd and even Kondo phases, which existed in the 
X = case. Qualitatively, the energetic preference of the 
odd Kondo phase can be seen from the approximate cor- 
respondence between the present model and the repulsive 
Hubbard model. Spin up (down) electrons in the Hub- 
bard model correspond to c (/) electrons here. In the 
Hubbard model, the dominant instability at half filling 
occurs at momentum (tt, tt) and corresponds to antifer- 
romagnetism. In the present model, antiferromagnetism 
clearly corresponds to the odd Kondo phase. 

In Fig. [3l we present the self-consistent numerical re- 
sults for the odd and even uniform Kondo phases. As 
one can see from Fig. [3 the quantum critical points sep- 
arating the FL* phase [6[ from the odd and even Kondo 
(heavy Fermi liquid) states occur at different values of 
interband coupling paramter V. Comparison of free en- 
ergies (see the inset of Fig. ^ , as expected, reveals that 
the odd Kondo phase has lower energy than the even 
phase. That the odd Kondo phase takes better advan- 
tage of the Fermi surface geometry is prominently man- 
ifested in the DOS. For the even Kondo state, although 
there is a depression exhibited in the projected conduc- 
tion electron DOS, the projected /-electron DOS shows 
a broad band behavior (see Fig. [IJa)). However, for the 
odd Kondo state, an s-wave-like gap opens in both the 
projected conduction electron and /-electron DOS (see 
Fig.Hb)). 

A finite chemical potential breaks the particle-hole 
symmetry in the conduction band and introduces a mis- 
match between the Fermi surfaces of the conduction elec- 
trons, FSc, and spinons, FS/. In this case, an instability 
to a state modulated at an incommensurate wavevector, 
related to the mismatch between FSc and FS/, is possi- 
ble. Indeed, in the continuum, an incommensurate har- 
monic Kondo wave has been recently found for circular 
Fermi surfaces jl2| . Here, we explore the real-space mod- 
ulation pattern, by solving the BdG equations ([5]) on a 
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FIG. 5: (Color) Spatial variation along the x-direction of the 
slave-boson b field (a), the effective hybridization hi{c[fi) (b), 
the averaged /-electron (c) and conduction electron (d) den- 
sities for V = 0.8 and ^l = -0.50. The temperature T = 0.01. 
The calculations are carried out on a lattice of 24 x 24 sites. 



finite size (24 x 24) lattice self-consistcntly via exact diag- 
onalization. To facilitate numerical convergence, we con- 
sidered only unidirectional patterns. This was achieved 
by "seeding" them as an initial guess, and iterating un- 
til convergence. We found that this procedure leads to 
convergence to a local minimum, without changing the 
initially imposed spatial period. The optimal modulation 
pattern is determined by comparing the free energies of 
the converged solutions with the real-space modulation 
periods N^jn^ (n^ is an integer running from to Nx ~ 1 
and Nx = 24 in our simulation). In general, the spa- 
tial modulation pattern is anharmonic, in close analogy 
to the stripes found in the mean-field treatment of the 
Hubbard model 



11| . For a fixed coupling parameter V, 



we find that there is a critical value of /i, below which a 
modulated Kondo state occurs. 

In Fig. [5] we show the spatial dependence of the slave- 
boson b field, the Kondo hybridization bi{c\fi), the lo- 
cal /-electron density, and the local c-electron density 
(the latter three quantities are gauge-invariant). We see 
clearly the emergence of an anharmonic Kondo stripe 
state when the system is deeply in the heavy fermion liq- 
uid phase. In this strong coupling limit, the hole density 
in the /-band is rather high, which leads to a short peri- 
odicity. We also notice that in the strong coupling limit, 
the "domain wall" between the regions of constant sign of 
hi is bond centered. We anticipate that by reducing the 
coupling strength V toward the Kondo breakdown crit- 
ical point, the anharmonicity will be suppressed along 
with a reduction of the modulation amplitude, leading 
to a nearly- harmonic Kondo wave [12l |. However, system 
size limitations and problems with numerical convergence 
in this parameter regime, prevented us from verifying 



this directly. The Kondo stripes are closely analogous 
to the antiferromagnetic (AF) stripes in the mean-field 
treatment of the repulsive Hubbard model [ll| . The field 
6(r) plays the role of the AF order parameter M(r), and 
the electron density n(r) corresponds to the hole den- 
sity p{v). The lowest-order symmetry-allowed coupling in 
the Kondo stripe free energy (which also follows from the 
constraint) is |6(r)|^n(r), which corresponds to S'^{y)p{y) 



term for AF stripes [21 1 . 

In summary, we have considered a lattice model of 
coupled conduction electrons and spinous, as a model 
of a heavy Fermi liquid (Kondo) state. We found that 
for nearest-neighbor hopping models, two distinct heavy 
Fermi liquid states are possible. We also point out that in 
the generic case of mismatched electron and spinon Fermi 
surfaces, electronically inhomogneneous heavy Fermi liq- 
uid phases (Kondo stripes) are expected to appear. 
These phases exhibit inhomogeneous Kondo hybridiza- 
tion and charge density, which should manifest in both 
local and bulk probes. Identifying these phases experi- 
mentally may help to determine the relevance of the spin- 
liquid physics to heavy fermion materials. 
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